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Abstract 

Contribution of nearby pulsars to the cosmic rays observed at Earth 
has been studied. It is found that the experimental bound on amplitude of 
cosmic ray anisotropy may produce significant constraint on the efficiency 
of converting pulsar rotational energy to emitted particles kinetic energy. 
Cosmic ray fluxes from two well known nearby gamma ray pulsars, namely 
the Vela and Geminga pulsars, are estimated. The analysis suggests that 
observed bound on cosmic ray anisotropy restricts the contributions of 
the Vela and the Geminga pulsars to at most 1 % of the observed cosmic 
rays below the knee. 
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1 Introduction 

Pulsars are widely considered to be natural sites for acceleration of charged 
particles (e.g. [1-12]). Along with the Supernova explosions and Gamma Ray 
bursts, they are regarded as the most probable sources of ultra high-energy 
cosmic rays within the galaxy. From the point of view of available energy, pulsars 
are even a more promising source because the rotational energy at birth can be 
more than ten times of the energy released in a typical supernova explosion. It 
also has the necessary conditions for accelerating particles to ultra high energies. 
Moreover observation of pulsed and/or steady flux of electromagnetic radiations 
from radio to gamma ray energies provides direct evidence that some pulsars are 
site of energetic particles of at least several TeV. Another important feature is 
the maximum energy (E max ) attainable by a particle in the acceleration process. 
In the case of SNR this is at most 10 15 eV [13]. On the other hand according 
to the Hillas condition [14], the maximum energy of a particle of charge Ze that 
can be contained near the light cylinder of a pulsar of angular speed Q rad s _1 , 
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radius R and with the surface magnetic field B s = B i2 x 10 12 Gauss is 

E max = 3.4 x l "Z* ia (^) 2 (^) 3 eV (1) 

which shows that the maximum possible energy can be very large, for fast 
rotating pulsars E max even could reach around 100 EeV [1,7], the highest energy 
cosmic ray particles observed so far. A point that usually has been raised 
against the pulsar model of origin of cosmic rays is that the derived spectrum of 
cosmic rays from pulsars is much flatter than the observed spectrum. However, 
such derivation is based on the assumption that the production spectra from 
all pulsars have the same slope and E max . It has been shown recently [9,10] 
that the expected cosmic ray spectrum in the energy range from PeV to EeV 
coincides with the observation if the distribution of pulsar initial periods is 
similar to the Gamma distribution [9] or if the log of initial pulsar periods and 
surface magnetic fields are given by the Gaussian distribution [10]. 
An important characteristics of cosmic rays is that they are highly isotropic; 
the amplitude of anisotropy is less then 10~ 2 below the knee energy region 
[15]. While estimating contribution of a pulsar such isotropic nature of cosmic 
rays has to be respected. But this point has not been taken into consideration 
in several investigations. In most of the recent efforts of estimating cosmic 
rays from pulsars collective contribution of all galactic pulsars are studied [7- 
10]. Hence it is not clear from those analysis whether the adopted cosmic ray 
production scenario in pulsars is consistent with the observed isotropy of cosmic 
rays or not. Here it is worthwhile to mention that the observed upper limit of 
cosmic ray anisotropy leads to some important conclusions on the supernova 
remnant (SNR) model of origin of cosmic rays as well as on their propagation. 
For instance in a recent work [16] it is shown that the observed cosmic ray 
isotropy does not support the idea the that only a small fraction of supernova 
remnants provide the main contribution to the cosmic ray flux in the knee region 
[17]. In another work [18] it is found that the SNR origin model of cosmic rays 
is compatible with the observed upper limit of cosmic ray anisotropy if the 
cosmic ray propagation is by diffusion with re-acceleration; the plain diffusion 
model of propagation gives too large anisotropy. Thus the study of individual 
contribution of pulsars, particularly from nearby pulsars, is of great significance. 
In the context of single source model of the knee [19] contribution of nearby 
pulsars to the cosmic ray spectrum at knee region has been studied recently [20] 
where parameters are chosen suitably so that the contribution of the pulsars 
becomes significant at the knee energy of the spectrum. But isotropic nature of 
cosmic rays was not incorporated in those studies. In the present work cosmic 
ray flux from nearby pulsars are estimated respecting the observed isotropy of 
cosmic rays. 

The plan of the work is as follows. After giving a short review of the acceleration 
of charged particles by pulsars, we obtain the cosmic ray spectrum of charged 
particles produced by a pulsar in sec. II. The contribution of a nearby pulsar to 
the cosmic rays at the Earth is estimated in section III. In section IV constraint 
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on the contribution of cosmic rays from a point source that may arise from 
the observed isotropy of cosmic rays are discussed. Fluxes of cosmic rays from 
two nearby gamma ray pulsars, the Geminga and Vela pulsars, are estimated in 
section V. Finally the results are concluded in section VI. 

2 Acceleration of nuclei by pulsars 

Pulsars arc generally believed to be rotating neutron stars [21]. The rotation 
of the magnetic dipolc of pulsars generates huge electrical differences between 
different parts of its surface and as a result the rotating neutron star should be 
surrounded with charged plasma, which is called a magnetosphcre. Since the 
moment of inertia of a neutron star is around 10 45 erg s 2 , a millisecond pulsar 
has a rotational energy E = \IQ 2 ~ 10 52 ergs. A fraction of such a huge rota- 
tional energy of a pulsar may be converted to the kinetic energy of the particles 
those present in the magnetosphere. 

The pulsar magnetosphere is usually considered to be composed of electron- 
positron pairs. They are expected to be produced to get rid of electric fields 
that are parallel to magnetic fields at a region where net charge density is dif- 
ferent from the Goldreich- Julian density [1]. However, hadronic component also 
may exist [1-12] in magnetosphere. Ions in the magnetosphere of a neutron star 
are most likely to come from the surface of the star. The surface of the neutron 
star is supposed to compose of heavy nuclei, most probably of iron nuclei [12]. 
The binding energy of these surface nuclei is, however, not accurately known. In 
the case of small binding energy (< 1 keV), thermionic emission of ions from the 
surface of the star is possible even for the old stars which may give practically 
free outflow of ions in the magnetosphere. But if the binding energy of surface 
atoms is large (2 — 3 keV), thermionic emission can take place only for young 
pulsars (surface temperature > 10 6o C). The local temperature of the polar cap 
can be higher due to the bombardment by particles of the magnetosphere and 
thus thermionic emission still may occur in pulsars of intermediate age. Surface 
atoms also can be stripped off by the strong electric field of the polar cap of a 
pulsar [4]. 

These nuclei can be accelerated by pulsars through large potential drop asso- 
ciated with strong electric field parallel to the pulsar magnetic field. Several 
detailed mechanisms have so far been suggested for accelerating particles by 
pulsars including the popular polar gap [4], and the outer gap [5] models. In 
the former model, acceleration of particles takes place in the open field line re- 
gion above the magnetic pole of the neutron star whereas in the case of outer 
gap model it occurs in the vacuum gaps between the neutral line and the last 
open line in the magnetosphere. Thus, the region of acceleration in the polar 
gap model is close to the pulsar surface, while the same in the outer gap model 
is close to the light cylinder. It has also been suggested that pulsars may accel- 
erate protons and heavy nuclei by converting their rotational energy to particle 
kinetic energy via a relativistic MHD winds near the light cylinder [3,7]. In 
this process particles may gain energy from the largest potential drop associ- 



3 



ated with a neutron star. In the present investigation, however, we restrict our 
discussion within somewhat general framework of electromagnetic acceleration 
process without focusing on any specific model. 



2.1 Maximum energy 

A large-scale electric field in the magnetosphere is expected to develop due to 
the rotating magnetic dipole of the pulsar. Neglecting the influence of the pulsar 
magnetosphere, the maximum electrostatic potential at time t is 

4>m = ^B s n 2 (t)R 3 (2) 

But the whole potential drop may not be available for the particle acceleration 
due to development of electron-positron pair cascades (originate from curvature 
radiation photons emitted by a seed electron accelerated by the drop) in opposite 
direction which quickly shuts off the potential drop [5, 12, 22]. 
If / fraction of the (j) m is responsible for the acceleration process, the energy of 
the accelerated nuclei would be 

Br-VTxtffEB*^)'^)' G.V. (3) 

In general / may depend on pulsar period and magnetic field. On the other 
hand if we consider the acceleration scenario in which magnetic field energy 
near the light cylinder is transferred to particle kinetic energy [3,7], the magnetic 

energy per ion at the light cylinder is E = 87m '^ - where Bi c = B s R 3 £l 3 /c 3 is 
the magnetic field at the light cylinder and riQj(r) = B s R 3 Vt/ (AirZecr 3 ) is the 
Goldreich- Julian density. Thus the numerical value of the magnetic energy per 
ion becomes the same as Eq.(3) with / = 1. 



2.2 Energy spectrum 

In the acceleration picture mentioned above it is natural to consider that at 
any evolutionary stage of the pulsar, all the particles are accelerated under the 
same potential difference and hence should attain the same maximum energy 
E cr . With the spin down of the pulsar, attainable energy of the accelerated 
particles changes. If we assume that the £ fraction of total rotational energy 
loss of a pulsar goes to accelerate nuclei then one can write the rate of number 
of accelerated particles at any moment 

However, the charged density follows from the above equation may not be always 
less than the Goldreich- Julian density (p w — ilB/(2nc)) [8]. But the expression 
for the GJ density is valid only in the co-rotation portion of the magnetosphere 
[1,8] whereas the most likely site of production of high-energy cosmic rays is 
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near the light cylinder. If r\ fraction of total rotational energy loss is due to the 
emission of magnetic dipole radiation, then 

So one can express the rate of change of rotational energy of a pulsar at any 
instant t as 

Ero< = 3.9 x lO^V 1 ^ (tS^tV (?£-)* GeV s" 1 (6) 



The evolution of the pulsar period is given by the expression 

PHU) = n tl) + 4.8 x 10- te - tl) B i2 (jp^-,) " ,7, 

Here P(t2) and P(ti) are the period of pulsar in seconds at time £2 and t\ 
respectively. The rate of change of pulsar period at a particular instant t is 
followed from the above equation which is given by 

With the evolution of the pulsar period, the energy of the accelerated particles 
changes according to Eq.(3) through Eq. (7). The resulting energy spectrum of 
the accelerated particles thus reads as 

dN 1.9xl0 45 / fr7n / R \ 3 \ 1 f rjl \ „ T . 1 

Apparently the energy spectrum of the produced accelerated particles has the 
E^ 1 behavior. But since the production times for particles of different energies 
are different, the observed spectrum may not go simply as E" 1 [8]. 



2.3 Escape from the nebula 

A young pulsar is generally encircled by the remnant of the pre supernova star. 
In fact most gamma ray pulsars have plerions as well as they have associated 
supernova remnants. Only the oldest pulsars such as Geminga or PSR1055-52 
show no sign of having plerions. So pulsar accelerated nuclei will inject into the 
nebula. Initially the accelerated nuclei will interact with the dense matter of 
the expanding envelope. The interactions, however, become negligible when the 
column density is very much less than the interaction length of the energetic 
nuclei in the medium. The column density of the matter in envelope offered 
to the injected high energy nuclei after t years from the birth of the pulsar is 
approximately 

YM M ™ b znf 1 Y 2 f M neb \ 2 f E SN V 1 _ 2 

() ^ =50 l^J ) U^J 9cm (10) 
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where M ne b and R ne b stand for mass and effective radius during the free ex- 
pansion phase of the envelope respectively. The radius of the envelope varies 
as R ne b — R + v t where R Q is the radius of pre-supernova star (R D ^ 10 14 
cm) and the velocity of the free expansion v ~ (2E SN /M neb ) 1/2 , the observa- 
tions suggest that the free expansion velocity of Crab nebula is around 2000 
KM s~ x . The condition for iron nuclei to pass through the envelope without 
major losses is that £ < 100 gcm~ 2 . Thus the energetic nuclei can escape the 
remnant without significant losses shortly after the explosion. 
However, energetic nuclei injected into the envelope may be trapped by the mag- 
netic field of the nebula. Outside the light cylinder the azimuthal component 
of the magnetic field dominates over the radial field and is given by 



where a is the ratio of the magnetic field flux to the particle energy flux at shock 
radius r. For the Crab pulsar a is estimated to be 0.003 [23] whereas for the 
Vela pulsar the value of a is close to 1 [24] . According to the standard picture 
of development of pulsar envelope [25] at the initial stage envelope expands 
freely with velocity v. At this stage the radius of the envelope increases with 
time roughly as 6.3 x 10 15 t cm where t is in years. With time the envelope 
slowly begins to decelerate by transferring its energy to the ambient medium. 
When the envelope reaches the Sedov phase, which occurs when the matter 
density of the nebula becomes comparable to that of the surrounding medium 
and begins at t sec iov — Rsedov/v where the radius of the envelope at Sedov 

1/2 

phase Rsedov — ( M ° n nv ) with n ~ 0.3 gcm~ 3 being the density of the 
surrounding medium (typically t se d ov ~ 2000 years), the velocity of envelope 
can be approximated as v se dov — v(t se dov/t) 3 ^ 5 [26]. Consequently the radius 
of the envelope goes with time as r ne b ~ v ^ledov^^ '■ During this stage the 
strength of the magnetic field in the outer part of the envelope takes a value 
around 3 to 5/iG, the typical value for the interstellar medium. Accordingly, the 
Larmor radius (rj, = E/{ZeB)) in the envelope becomes about 10 18 ( 1 p eV ) ^ 
cm. Thus, nuclei of energy even 1 PeV should be trapped by the magnetic 
field of the nebula. However, it is known that several instabilities develop in 
the process of confinement of charged particles by magnetic fields and the high- 
energy nuclei will finally escape along the field lines of the irregular fields [8] . 
The energetic nuclei will propagate diffusively in the envelope before escaping 
into the interstellar medium provided diffusion time Tdiff (= r 2 /AD(E), D(E) 
is the diffusion coefficient) for traversing the radial distance r„ e & is much larger 
than t = r ne b/c, the time required for straight line propagation. When the 
energy of particles is too high they move on almost rectilinear trajectories. The 
pulsar accelerated particles will escape from the nebula when the mean radial 
distance traveled by the particles becomes comparable with the radius of the 
nebula at the time of escaping. The mean radial distance traveled in diffusive 




(11) 
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propagation in time t is [27] 



< r rad >= 2 ] J^ (12) 

where D{E) is assumed as time independent. Such an assumption holds, at 
least as first approximation, only after the nebula enters the Sedov phase, as 
thereafter the magnetic field of the nebula does not change appreciably with 
time. The Diffusion coefficient for relativistic particles is defined by the relation 
D(E) = icAd, where Ad is the mean scattering length by magnetic irregularities. 
The problem is that Ad is not known with good accuracy. Theoretical limit on 
Ad is that it cannot be smaller than the Larmor radius [28] . Thus the minimum 
value of diffusion coefficient (so-called Bohm diffusion constant) corresponding 
to maximally turbulent medium is given by 

D = X -cr L . (13) 

But after the nebula enters into the Sedov phase the magnetic field fluctuations 
is more likely to represent by the Kolmogorov spectrum [12]. In that case the 
diffusion coefficient is well described by the relation 

D(E) = \rT?Jl ■ (14) 

where l ce ii represents the characteristic size of the magnetic cells (irregularities) 
in the nebula which is obviously less than the radius of the nebula but can be 
within an order of that, the characteristic size of filamentary structure of the 
Vela supernova remnant as revealed from multiwave length studies indicates 
such possibility. 

It is evident from the foregoing discussion that mean storage time of cosmic rays 
of different energies in the nebula depends on the choice of the diffusion coef- 
ficient. But even for the Bohm diffusion, which gives maximum confinement, 
mean storage time of a PeV proton in the nebula is only few hundred years 
whereas the same for a PeV iron nucleus is few (typically 6 to 10) thousand 
years. 



3 Cosmic ray flux from a pulsar to be observed 
at Earth 

The diffusion process governs the propagation of accelerated charged nuclei from 
the source. Usually the diffusion scenario is considered as Gaussian. Neglecting 
the effect of energy gained/losses during propagation, convection, losses of nuclei 
by collision and nuclear interactions, the diffusion equation is given by 

— =A.(D(E)AN(r,t,E)+Q(r,t,E) (15) 
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where Q(r, t, E) represents the source term. The Green's function for the above 
equation is given by 

° = 8(nDt)^ eXp[ - r2/{ADt)] (16) 

Here t is the time when cosmic ray of energy E escapes from the pulsar nebula. 
Thus the intensity of cosmic rays of energy £ at a distance r from a pulsar 
would be 

I cr (r,E) = £2.6 x 10- [fZB 12 m {-^)) ( W ^ ? ) G(r,r) 

cm^s^sr^GeV- 1 (17) 

where c is the speed of light. 



4 Anisotropy 

The anisotropy predicted for galactic cosmic ray sources relies on the model for 
the production of cosmic rays and for their propagation. The production of 
cosmic rays depends on the nature of the sources whereas the diffusion process 
due to scattering from minor irregularities in the field governs the propagation 
of accelerated charged nuclei from the source to Earth. The resultant amplitude 
of anisotropy is given by [29] 

s- Xd \ dI{E) \ ri8i 

where I is local cosmic ray intensity and dl/dr is the intensity gradient. Em- 
ploying the simplified model of diffusion as mentioned in previous section, one 
gets 

'-^ssfe (19 » 

Here h(E) denotes the ratio of the cosmic rays of energy E from the source 
to the total observed flux of cosmic rays at the same energy from all sources. 
The energy dependence of S relies mainly on h(E). But emission time and 
hence t also can be different for cosmic rays of different energies depending 
on the nature of source. A nice feature of the expression (19) is that it does 
not depend on the diffusion coefficient. However, if a different propagation 
scenario, such as the so-called anomalous diffusion [30], is adopted then the 
expression for anisotropy amplitude might depend on the diffusion coefficient. 
Here it should be mentioned that the numerical value of diffusion coefficient 
at concerned energy range is quite uncertain. Another important point is that 
once the contribution of the source to the total cosmic ray flux is fixed, there is 
no adjustable parameter left in the expression (15). 

As the observations show that cosmic rays are highly isotropic, S is restricted to 
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a small value. For example harmonic analysis of the right ascension distribution 
of cosmic rays gives amplitude of anisotropy upto the knee region is around 10~ 2 
to 10~ 3 [15]. The measurements at the higher energy region suffer from the low 
flux of cosmic rays and thus S is less restricted there. While estimating flux 
from a point source constraint imposes by the observed anisotropy through the 
Eq.(19) has to be respected which in turn may restrict the parameter £. The 
phase of the amplitude of first harmonic is also a quantity of significance; the 
contribution of a point source to the total cosmic ray intensity needs to be 
consistent with the observed phase of the amplitude [15]. 

5 Application: Cosmic rays from two nearby 
gamma ray pulsars 

In a recent work [8] it has been shown that Gamma ray pulsars of the Galaxy 
could contribute a significant fraction of the total cosmic ray flux at and above 
knee energy region. However, the contribution of the nearby sources has not 
been considered there. Effect of observed bound of cosmic ray anisotropy on the 
parameters involved in acceleration of cosmic rays by pulsars, if any, is expected 
to reveal from the study of contribution of cosmic rays from the nearby pulsars. 
For pulsars of large distances the effect of different emission (form pulsars) 
times of cosmic ray particles having different energies on the observed cosmic 
ray spectrum are not significant. This is because in such cases the propagation 
time is large compared to the life time of the pulsars. The effect of different 
time of emission could be important only for nearby pulsars. In the following 
cosmic ray fluxes from the two nearby gamma ray pulsars, namely the Vela and 
the Geminga, are considered. The present analysis is restricted to gamma ray 
pulsars only as the emission of gamma rays is a definite signature of presence of 
high-energy particles. We take D(E) = D (E/Z) S cm 2 s~ 1 (where E is in GeV) 
with D = 1.4 x 10 28 for galactic magnetic field and 5 = 1/3 [31] . 

5.1 The Vela Pulsar 

Vela was found by SAS 2 as the brightest object in the gamma ray sky. It is 
comparatively close to the Earth and so the surrounding nebula is well studied. 
Its characteristic age is around 10, 000 years with periodicity 89.3 ms, slow down 
rate is 1.25 x 10~ 13 , and surface magnetic field is around 3.4 x 10 12 Gauss [32]. 
The distance of the object from the Earth is around 500 pc [33] though recent 
works suggest a smaller value of 300 pc [34]. 

The observed pulsar period and slow down rate are consistent with both the 

Eqs. (5), and (8) for (jq^) 6 ( 10 4b"/ cm . ) = 0.24. The initial pulsar period 
is estimated as 26.3 ms from Eq.(7). In the framework of the self-sustained 
outer gap model [35] in which / is given by 

/^l.OTxlO-^^)" 2721 ^ 877 , (20) 
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Figure 1 : Estimated cosmic ray flux from the Vela pulsar 



Figure 1: Estimated Cosmic Ray flux from the Vela pulsar 



the maximum energy would be 30 Z TeV. Since the pulsar is relatively young 
in the case of Bohm diffusion it is likely that the pulsar accelerated cosmic rays 
still remain confined in the nebula. But if the characteristic size of the magnetic 
irregularities in the nebula is larger than the Larmor radius then the diffusion 
coefficient will be given by the Eq.(14) and assuming l ce u is of the same order to 
the radius of the nebula (corresponds to minimum confinement) it is found that 
the mean storage time of a 30 TeV proton is around 3000 years; the same for a 
100 TeV iron nuclei is about 6000 years. Resulting spectra of cosmic rays from 
the pulsar to be observed at earth are shown in figure 1 for both proton and 
iron primaries assuming £ = 0.01 and 0.001 respectively. The parameter £ is so 
chosen to respect the observed isotropy of cosmic rays. The observed all particle 
cosmic ray (best fitted) spectrum is also given in the figure for comparison. 
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5.2 The Geminga pulsar 



Another pulsar in the area is Geminga (PSR J0633+1746), which is discovered 
by the SAS-2 group [36] and later confirmed by the COS-B group [37], about 
150 pc away from the Earth. Its radial velocity is unknown, but if it were 
200 KM/sec, it could have been within 100 pc of Earth at 340,000 years ago. 
Geminga is a unique object: a highly compressed, spinning neutron star which 
does not emit radio waves like the other well-known pulsars. Yet it is a powerful 
source of pulsating gamma-rays and X-rays. Geminga is now known to be a 
rotation-powered pulsar with period P = 0.237 s, P = 1.0975 x 10~ 14 and 
surface magnetic field B = 1.6 x 10 12 G [38]. 




Figure 2: Estimated Cosmic Ray flux from the Geminga pulsar 

Since Geminga is a relatively old pulsar, presence of ions in its magneto- 
sphere is certainly questionable. However, observation of high-energy gamma 
rays from this pulsar indicates the presence of energetic particles in its atmo- 
sphere. Besides the non-observation of radio waves from the pulsar suggests 
that these high-energy particles may not be electrons. Hence it is likely that 
ions are present in its magnetosphere; at least one cannot rule out this possibil- 
ity. Even if sufficient ions are not available in the pulsar atmosphere presently, 
one needs to consider the high energy particles which were produced at earlier 
stages. According to Eq.(7), the initial pulsar period is 36.6 ms. Here also we 

take ( 10 e fl cm ) 6 ( w^gem? ) = °- 24 wmcl1 follow from Eq. (5) and is consis- 
tent with Eq. (8). The maximum energy of the energetic nuclei thus would be 
only 4Z x 10 13 eV when / is given by Eq.(20). Taking £ = 0.1 so that observed 
isotropy is mostly obeyed, we estimate the spectra of cosmic rays from the pul- 
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sar to be observed at earth for proton and iron primaries which are shown in 
figure 2. Observed (best fitted) all particle spectrum [39] is also presented for 
comparison. It is worthwhile to mention that for Geminga the nature of dif- 
fusion process of cosmic rays in the supernova envelope has little effect on the 
resultant (observed) spectra because of its age. 

6 Discussion 

In the present work contribution of a nearby pulsar to cosmic rays observed 
at earth has been studied respecting the observed isotropy of cosmic rays. It 
is found that the parameters involved in the cosmic ray generation process 
in pulsars, particularly the efficiency of converting pulsar rotational energy to 
emitted particles kinetic energy may be severely constrained by the observed 
limit on cosmic ray anisotropy. Consequently production of high-energy gamma 
rays through collision of cosmic rays with the ambient matter would be also 
restricted [40]. 

Cosmic ray fluxes from two well-known gamma ray pulsars, namely the Vela and 
Geminga pulsars, are obtained. It is found that owing to the observed bound 
on cosmic ray anisotropy the contributions of the Vela and the Geminga pulsars 
are limited to at most 1 % of the total cosmic ray intensity below the knee. This 
means the Vela and the Geminga pulsars do not contribute substantially to the 
total cosmic ray intensity. Such a conclusion is consistent with the observed 
phase of the amplitude of first harmonic; below the knee energy maximum flux 
of cosmic rays (phase of the anisotropy) is observed from the second quadrant 
of the Galaxy [15], while both Vela and Geminga are in the third quadrant. 
Another interesting observation is that the cosmic ray spectra to be observed at 
Earth from both the pulsars, the Geminga and the Vela, are found noticeably 
different than the production spectra which is mainly due to combined effect 
of different generation (and emission) time of pulsar accelerated cosmic rays of 
different energies and their propagation. 
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